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ABSTRACT 

One of the goals in quantum simulation is to adiabatically generate the ground state of a 
complicated Hamiltonian by starting with the ground state of a simple Hamiltonian and slowly 
evolving the system to the complicated one. If the evolution is adiabatic and the initial and final 
ground states are connected due to having the same symmetry, then the simulation will be 
successful. But in most experiments, adiabatic simulation is not possible because it would take 
too long, and the system has some level of diabatic excitation. In this work, we quantify the 
extent of the diabatic excitation even if we do not know a priori what the complicated ground state 
is. Since many quantum simulator platforms, like trapped ions, can measure the probabilities 
to be in a product state, we describe techniques that can employ these measurements to 
estimate the probability of being in the ground state of the system after the diabatic evolution. 
These techniques do not require one to know any properties about the Hamiltonian itself, nor to 
calculate its eigenstate properties. All the information is derived by analyzing the product-state 
measurements as functions of time. 

Keywords: Quantum simulation, ion trap, adiabatic state preparation, transverse field Ising model, ground state probability 


1 INTRODUCTION 

The Hilbert space that describes a strongly correlated many-body quantum system grows exponentially 
in the number of particles N, so determining the ground state of a complex many-body quantum system 
becomes numerically intractable when the size of the quantum system becomes too large to be represented 
on a classical computer (unless there is some other simplification, like weak entanglement, etc.). The idea 
to simulate complex many-body quantum systems on a quantum computer was proposed by Feynman in the 
i980’s m. Since Feynman’s proposal quantum alogorithms have been proposed to calculate eigenvalues 
and eigenvectors of these intractable systems Q.One of the challenges with creating the ground state on a 
quantum computer, say by adiabatic evolution of the ground state from a simple to a complex one, is how do 
we determine the extent of the ground-state preparation. After all, we don’t know what the ground state is a 
priori so it is difficult to know what the probability to be in the ground state is. In this work, we propose one 
method to determine the probability to remain in the ground state. While this analysis is applied to ion-trap 
emulators (that model interacting spin systems), the general discussion can be applied to any quantum 
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computer that performs ground state preparation, but creates diabatic excitations as a result of a too rapid 
time evolution. To date, trapped-ion quantum simulators have seen success in two different platforms: 
the Penning trap has trapped 300 ions in a planar geometry and generated Ising spin interactions ^ 
and the linear Paul trap has performed quantum simulations with 18 ions in a one-dimensional crystal 
The suecess of these traps as quantum simulators is attributed to their long coherence times, precise spin 
control, and high fidelity. Here, we will focus on the linear Paul trap quantum simulator. 

In an ion trap quantum simulator, hyperfine states of the trapped ions are used for the different spin states 
(for simplicity, we can consider only two states, and hence, a spin-one-half system). Optical pumping can 
be employed to create a product state with the ions all in one of the two hyperfine states with fidelities 
close to 100%. A coherent rotation of that state can then be used to create a wide range of different initial 
product states. By turning on a large magnetic field, this state can configured to be the ground state of the 
system. Then the magnetic field is reduced slowly enough for the system to remain in the ground state 
until the system evolves into the complex spin Hamiltonian in zero magnetic field. The challenge is that 
the evolution of the system must be completed within the coherence time of the spins, which often is 
too short to be able to maintain adiabaticity throughout the evolution (and indeed, becomes increasingly 
difficult as the system size gets larger). We propose to employ the time evolution of an observable, 
0{t) = {xp{t)\0\xp{t)) (with I'lpit)) the quantum state of the system at time t > tgtop)^ after the evolution 
to the final Hamiltonian, to measure the absolute probability of the ground state. The time evolution of the 
observable, 0{t), oscillates at frequencies given by the energy differences between the final eigenstates 
(where the Hamiltonian becomes time independent). More concretely, let H{tstop)\'m) = Em\m), then the 
time-dependent expeetation value satisfies 

0{t > tstop) = P^Pn{m\d\n) exp[-i{En - Em)t], (1) 

mn 

where Pm = {rn\'ip[tstop)) is the overlap of the state \'ip[tstop)) with the eigenstate \m) (we have set 
h = 1). Previously, we showed how Fourier transforming the time series and employing signal proeessing 
methods like compressive sensing, allows one to extract the energy differences as a type of many-body 
eigenstate spectroscopy iSlIll. Here, we focus on the amplitude of the oscillations, given by PiPn, which 
is proportional to the probability amplitude of the ground state (Pi), if the ground state still has a high 
probability in \^p[tstop)) - Note that we do not need to know the explicit ground-state wavefunction to extract 
its probability from these oscillations. This is the main advantage of this technique. 

We illustrate how the the ground state probability can be extracted by analyzing the amplitude of 
the oscillations of the simplest time-dependent Hamiltonian: the two-level Landau-Zener problem. The 
Landau-Zener problem is defined via 


n[t) = B^[t)a^ + aP (2) 

Here cr“ are the Pauli spin operators in the a = x, y, or z direction. The Pauli spin operators have the 
commutation relation 

where the Greek letters represent the spatial directions and is the antisymmetric tensor. The Landau- 
Zener problem has a minimum energy gap occurring when B^[t) = 0, as shown in Fig. Since the 
Landau-Zener problem is a two state system, the probabilities. Pi and P2, are related by P^ = 1 — Pf and 
the state, \'ip[tt)), can be represented by Pi = cos(0) and P 2 = sin(0). Using this state in Eq. ([^, we find 




(3) 
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Figure 1. (Color online.) The energy spectra of the Landau-Zener problem. There is a minimum energy 
gap when the excited state, in red, lies at its minimum energy difference with the ground state. This occurs 
at = 0. 

that the expecation value 0{t> tstop) becomes (neglecting terms with no time dependence) 

> tstop) = cos(0)sin(0) |(1|0|2} exp[-i (£’2 - Ei)t] + {2\d\l) exp[-i{Ei - £’ 2 )t]| • (4) 

The ground-state probability [cos^((/))] can then be calculated if (1|(!)|2) is known. The amplitude of the 
oscillations is sin(20)(l|(2|2). However, even if 0 is extracted from the amplitude, there are two always 
solutions, except when 0 = 7r/4 (see Fig.|^, and hence two possible ground state probabilities. In Fig.|^ 
we show both the amplitude of the oscillations and the ground-state probability as a function of 0, where 
the dashed line shows that the amplitude is not unique to a single ground-state probability. However, once 
the system has more states, one does not have a simple closed set of equations and the analysis of the 
amplitude of the oscillations can only deduce the ground-state probability when the ground-state amplitude 
is dominant in \tp{t)). We demonstate this below with the transverse field Ising model. 

It is well known that the amount of diabatic excitation in the Landau-Zener problem increases the faster 
the magnetic field is ramped from — oo to -foo. The general protocol that we employ is as follows (and is 
depicted schematically in Fig.|^: 

1. Initialize the state in an arbitrary state, in the following examples we initialize the state in the ground 
state of the Hamiltonian with a large polarizing magnetic field. 

2. Decrease the magnetic field as a function of time to evolve the quantum state, as shown in Fig. |^A), 
where, for concreteness, we show an example of a magnetic field that changes linearly. 

3. Hold the magnetic field at its final value which is first reached at t = tstop until the measurement is 
performed at the time interval tmeas. after the field has been held constant [see Fig.|^A)]. 

4. Measure an observable of interest, 0{t), for a number of different tmeas. values. 

5. Determine the amplitude of the oscillations. 

Note that one requirement of this approach is that the observable of interest must oscillate as a function in 
time as given in Eq. ([T]). The amplitude is extracted from the first maximum and minimum of the observable 
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Figure 2. (Color online.) Analytic solutions to the Landau-Zener problem given ground-state probability 

(black solid line) and the amplitude of the oscillations, P 1 P 2 , of 0{t > tgtop) (red solid line). The 
amplitude of the oscillations is not unique to a single ground-state probability, as highlighted by the dashed 
line. 


as a function of time by 


Amplitude = 


max[0(f)] — min[C>(t)] 


(5) 


The time evolution of the wavefunction \tp{t)) is calculated by solving the time-dependent Schrodinger 
equation: 

8 

( 6 ) 


by using the Crank-Nicolson method to time evolve the state IV'(^))- This technique solves the problem 
with the following approach |I7]] 


I + i^nit + 5t)j + 5t)) = (^- 

Note that the Hamiltonian is time-dependent until tgtop is reached, when it becomes constant in time. 


(7) 


We present a numerical example to illustrate this protocol by analyzing the oscillations for the Landau- 
Zener problem. Due to the fact that the eigenstates for the Landau-Zener problem at \B^{tstop)\ 3> 1 
approach the eigenstates of the cr* operator, if one measures an operator that is diagonal in this basis, there 
will be no oscillations in the expectation value. Hence, we measure the expectation value of the operator 
0(0), the Pauli spin matrix that points in the 0 direction. 


d{e) = R^{e)a^R{e), 

where R{d) is the global rotation about the y-axis and is given by 

R{0) = i cos f ^ j + ^ ^ ’ 


( 8 ) 


(9) 
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Figure 3. (Color online.) Schematic diagram of the experimental protocol. (A) The transverse magnetic 
field as a function of time is diabatically ramped down to a chosen value, B^itstop)- B^{tstop) is then held 
for a time interval of tmeas.- (B) An observable is measured at the end of the interval tmeas.- The amplitude 
of the oscillation is taken, for simplicity, to be the amplitude of the initial oscillation after tgtop- 


where 9 = 7 r /2 produces 0{9 = 7 r/ 2 ) = d®. 

For our numerical examples with the Landau-Zener problem, we use a linear ramp, B^{t) = rt + Bq, 
where 5o < 0. |5o| is chosen to be large in comparison to 1 to polarize the spin. We evolve the state to 
tstop, such that B^itstop) ^ 1- We present the time evolution for 4 different r = 1.5, 3.25, 5.0, and 9.0 for 
3 different 6* = 7r/9,7r/3, and 7r/2 in Fig.|^ The amplitude of the oscillations becomes 1 when 6* = 7r/2. 
When T = 5.0 the amplitude of the oscillations is maximized in comparison to the 3 other r’s. 

In Fig. 1^ we show the probability of the ground state compared to the amplitude of the oscillations as 
a function of r. The amplitude of the oscillations increases as 6 is increased. This is due to the term 
dominating the 0{9) operator instead of the term. When the ground state probability approaches 0.5 
the amplitude of the oscillations is maximal and it decreases either when the ground state probability 
increases or decreases. As the probability of ground state approaches 1 the amplitude is expected to become 
0, which can be obscured due to experimental noise. In order to determine which side of the maximum 
the measurement of the amplitude is on, multiple r experiments must be run to track the depletion of the 
ground state as the amplitude of the oscillations reaches a maximum and then decreases. Note further 
that in this case, since there is only one excited state, one can, in principle, always determine the ground- 
state probability by measuring the amplitude of the oscillation and extracting the appropriate probability 
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Figure 4. (Color online.) Time evolution of the measurement of 0(9) for r = 1.5 (black), 3.25 (red), 5.0 
(blue), and 9.0 (green) using 9 = irjQ (A), 7r/3 (B), and 7r/2 (C). The measurement is being done when 
B^{t) 3> 1. The amplitude of the oscillations converge to 1 as 6* is increased to 7r/2. 


amplitude. For more complex systems, such a procedure will not be possible, but the monotonic nature of 
the curve (at least while the probability for the ground state remains above 50%) will allow us to determine 
whether a given run of the experiment increases the probability to be in the ground state, which can be 
employed to optimize the ground-state preparation if it is done with some alternative quantum control 
method besides adiabatically evolving the system. Indeed, we believe this has the potential to be the most 
important application of this approach. 



Figure 5. (Color online.) Analysis of the amplitude of the oscillation as a function of the r for 6 different 
9"s. The dashed line is when the ground-state probability is 0.5. As 9 increases to 7r/2 the amplitude 
scale increases as well. When the ground-state probability is near 0.5 the amplitude of the oscillations is 
maximized and the amplitude decreases when either the ground state probability decreases or increases. 
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This simple example shows us a number of important points. First, one may need to rotate the 
measurement basis if the final product basis are eigenvectors of the Hamiltonian. Second, as the ground 
state is depleted, the amplitude of the oscillations grows until it reaches a maximum, when the system is 
equally populating both eigenstates. If the ground state is further depleted, the amplitude of the oscillations 
will decrease. One can make a mistake in estimating the probability in the ground state if one does not 
know which side of the curve one is on (probability of the ground state below or above 50%). On the other 
hand, if one knows which side of the curve one is on, due to making measurements at earlier times to track 
the ground-state depletion, then one might be able to further use the amplitude to determine the ground 
state probability in the Landau-Zener problem. When we change to the ion-trap system and examine the 
transverse-field Ising model, then the procedure becomes more complicated because there are more states 
that the ground state can be depleted into, and this complicates the analysis. 


2 TRANSVERSE FIELD ISING MODEL 

Now we describe a more realistic case of the transverse field Ising model. The transverse field Ising model 
for N particles is given by 


N N 

Hit) = -j±J2 E 

i<j 


( 10 ) 


where the Jy are the spin-spin interactions produced by a spin-dependent force and given by |l8]] 


N 




j/=i 


O O ) 


( 11 ) 


where hi^ is the normalized eigenvector of the phonon modes, is the corresponding frequency 
of the phonon mode, Vt is the single spin flip Rabi frequency, and vji is the recoil frequency associated 
with the dipole force on an ion, from which we define our energy units with Jq = . The Raman 

beatnote frequency /i is tuned to the blue of the largest ujp (which here is the center-of-mass phonon, 
cocom)- The details of calculating the Jij can be found elsewhere |l9l|. The oc \rij\~°‘ with the 
interparticle distance and the exponent a being tunable between 0 and 3. The exponent a is tuned by 
changing /r or by changing the ratio of the longitudinal to the transverse trap frequencies. Here, we study 
the ferromagnetic interaction of the Ising model with J± > 0. The Pauli spin matrices are now associated 
with each lattice site. The of the transverse Ising model have spatial-reflection symmetry such that 
Jij = Jji and the eigenstates have the same symmetry. The eigenstates of the transverse field Ising model 
also have spin-reflection parity, that is, under the partial inversion transformation cr* —?• cr*, 
and cr^ —> — cr^. The spin-reflection parity and spatial-reflection symmetry produce avoided crossings 
between eigenstates with the same parity and symmetry, such that a minimum energy gap to the lowest 
coupled state occurs as shown in Fig.|7j 


The experimental protocol is essentially the same as before, except for a few differences. The first 
difference being that the transverse magnetic field now depends exponentially as a function of time and is 
given by 


B^{t) = Boew *, 


( 12 ) 
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Figure 6. (Color online.) Schematic diagram of the experimental protocol used for the transverse field 
Ising model. (A) The transverse magnetic field as a function of time is diabatically ramped down to a chosen 
value, B^{tstop = 6r). Then B^{t) is quenched to 0 before the start of the measurement time interval, 
tmeas.- (B) The transverse magnetic field is ramped down and held constant at a final value determined by 

t — tstop- 


as shown in Fig.|^A). The second difference is that we will perform two different experimental protocols 
where the initial state is evolved to tgtop = 6^ and before the time interval tmeas. starts, the magnetic field 
is quenched to zero B^(tmeas.) = 0, as shown in Fig.|^A), or the transverse magnetic field is held at its 
final value which is first reached at t = tstop, as depicted in Fig. |^B). We work with parameters for the ion 
chain where the exponent CK 1. The energy spectra are plotted in Fig. |7| 



Figure 7. (Color online.) The energy spectra of the ferromagnetic transverse field Ising model with N=10. 
The red curve shows the first coupled excited state and the minimal gap occurs at B^ 0.72. 
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There are a number of additional complications. First off, the eigenstates at S® = 0 are product states 
along the 2 ; direction, hence we need to rotate again to see the oscillations. We choose 0(9) to be the 
average magnetization in the 6*-direction 


= (13) 

i—1 

where R{d) is now the global rotation given by 

N 

R(D )=n 

i—1 

where 9 = 2 yields crf^^. Measuring the average magnetization in the 6*-direction produces the needed 

oscillations. 


, ^ \ ^ 

I cos - + iaf sin - 

2/ * V2 


(14) 



Figure 8. (Color online.) Three different examples of 0{9) as a function of time when r = 0.2 (black), 
0.4 (red), and 0.6 (blue) for 9 = tt/6 (A), tt/S (B), and ir 12 (C) at B/Jg = 0. 

We next show simulated data for the transverse field Ising model with J± = 1 and Jq = IkHz. The 
parameters for the are // = 1.0219ucom and the antisymmetric ratio of the trap frequencies is 0.691/4.8 
which results in an a 1.0. The initial state is evolved to tgtop = 6 t and before the time interval tmeas. 
the field is quenched to zero, as shown in Fig. |^A). 

In Fig.[^ we show the time evolution of 0{9) with 9 = vr/O, tt/S, and 7r/2 for 3 different tJq = 0.2, 
0.4, and 0.6. The amplitude of the oscillations follow a similar trend to Fig. Q such that at r = 0.4 the 
amplitude of the oscillations are at a maximum in comparison to other r’s. Additionally, as 9 is increased 
to 7r/2 the amplitude of the oscillations increase as previously seen in the Landau-Zener example. 

In Fig.|^ we compare the probability of the ground state to the amplitude as a function of the ramping r. 
In general, the amplitude of the oscillations is maximized near r = 0.4 when the ground state probability 
is ~ 0.61 and the amplitude decreases as the probability to be in the ground state either increases or 
decreases. Similar to the Landau-Zener problem, as the probability to be in the ground state increases 
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to 1 the amplitude will decrease to 0. However when the ground state probability approaches 0.5, the 
amplitude of the oscillations is at a local minimum. As previously seen in the Landau-Zener example, a 
single measurement of the amplitude ramped at r cannot determine whether the probability of the ground 
state is high or low. Hence, the probability needs to be tracked by using a series of measurements. 



Figure 9. (Color online.) Amplitude of the oscillations at S/Jo = 0 as a function of different r as 
compared to the ground state probability. The arrows for each curve point toward the appropriate vertical 
axis. As the 9 increases, the amplitude of the oscillations increases as well. The amplitude of the oscillations 
become a maximum when the ground state probability is near 0.6 and the amplitude decreases when the 
ground state probability increases above or decreases below 0.6. When r = 0.3 there is a local minimum 
in the amplitude. 


In general, the analysis of the amplitude can be done for different tgtop’s in which the transverse magnetic 
field is held constant at the strength of B^{tstop), shown in Fig.|^B). Fig. 10 shows the amplitude of the 
oscillations as a function of B^{tstop) for 3 different r = 0.2, 0.4, and 0.6. The amplitude of the oscillations 
increases as the transverse magnetic field approaches the minimum energy gap and the excitations are 
created from the ground state, depending on the r. Once past the minimum energy gap, the ground-state 
probability increases as de-excitations occur and conversely the amplitude of the oscillations decrease as 
well. However a similar response will occur if excitations are being created after the minimum energy gap. 
Unfortunately, the analysis of the amplitude will not distinguish between these two possibilities. 


3 CONCLUSION 

In this work, we have proposed to analyze the amplitude of the oscillations for a given time-dependent 
Hamiltonian that is held constant for a time interval tmeas. to extract information about the ground-state 
probability. We demonstrated this analysis for the Landau-Zener problem and for the transverse field 
Ising model (as would be simulated in the linear Paul trap). In both of the Hamiltonians, the amplitude 
of the oscillations becomes a maximum at a particular probability of the ground state and decreases as 
the ground state probability either increases or decreases. Hence a single measurement of the amplitude 
cannot determine which side of the maximum one is on. Therefore multiple measurements are needed to 
be made where the amount of excitations are varied. Additionally as the probability of the ground state is 
approaching 1, the amplitude decreases to 0 which can be difficult to measure given experimental noise. 
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Figure 10. (Color online.) Amplitude of the oscillations when the transverse field Ising model is held at 
different 5/ Jo for 3 different r = 0.2 (black), 0.4 (red), and 0.6 (blue), (where the arrows point toward the 
appropriate vertical axis). Before the minimum energy gap, the oscillations increase as the probability of 
the ground state decreases. However, after the minimum energy gap is passed, depletion of the excited 
states back to the ground state occurs and the amplitude of the oscillations decrease accordingly. This 
depletion is difficult to detect by measuring the amplitude of the oscillations at only one time. 


In this work, we have described the simplest analysis one can do to extract information about the 
probability of the ground state. This approach can be refined by using signal processing techniques like 
compressive sensing ifT^ to determine the Fourier spectra of the excitations. By monitoring the change 
of the weights of the delta functions, one can produce more accurate quantitative predictions for the 
probability of the ground state, because we can directly measure Pf Pm for a few different m values. But 
this goes beyond the analysis we have done here. 

For the transverse field Ising model, de-excitations are observed, and were reflected in the amplitude of 
oscillations. However, after the minimum energy gap, more diabatic excitation can be created, but it is 
difficult to distinguish between the de-excitations and excitations. One interesting aspect is that as long 
as the ground-state probability remains high enough, measuring the height of the oscillation amplitude 
can be used to optimize the ground-state probability as a function of parameters used to determine the 
time-evolution of the system. This can be a valuable tool for optimizing the adiabatic state preparation 
protocol over some set of opimization parameters. 
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